Creep rupture of viscoelastic fiber bundles 
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We study the creep rupture of bundles of viscoelastic fibers occurring under uniaxial constant 
tensile loading. A novel fiber bundle model is introduced which combines the viscoelastic constitutive 
behaviour and the strain controlled breaking of fibers. Analytical and numerical calculations showed 
that above a critical external load the deformation of the system monotonically increases in time 
resulting in global failure at a finite time tf, while below the critical load the deformation tends to a 
constant value giving rise to an infinite lifetime. Our studies revealed that the nature of the transition 
between the two regimes, i.e the behaviour of tf at the critical load o c , strongly depends on the 
range of load sharing: for global load sharing tf has a power law divergence at <r c with a universal 
exponent of 0.5, however, for local load sharing the transition becomes abrupt: at the critical load 
tf jumps to a finite value, analogous to second and first order phase transitions, respectively. The 
acoustic response of the bundle during creep is also studied. 
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Fiber reinforced composites are of great technologi- 
cal importance due to their very good performance un- 
der extreme circumstances. Under high steady stresses 
these fiber composites may exhibit time dependent fail- 
ure called creep rupture, which limits their life time and 
consequently has a high impact on the applicability of 
these materials for construction elements. Both natu- 
ral fiber composites like wood EU2] and various types 
of fiber reinforced composites |3-|5| show creep rupture 
phenomena, which have attracted continuous theoreti- 
cal and experimental interest over the past years. The 
underlying microscopic failure mechanism of creep rup- 
ture is very complex depending on several characteristics 
of the specific types of materials, and is far from be- 
ing well understood. From technological point of view, 
one of the most important aspects of creep rupture is 
the statistics of life time (or time to failure) as a func- 
tion of the external steady load, however, only a limited 
number of systematic experimental works is available for 
fiber reinforced composites q-h], more information has 
been accumulated about natural fiber composites 
In Ref. § a theoretical model of creep rupture of brittle 
matrix composites reinforced with time dependent fibers 
was worked out, where the fibers are assumed to have a 
finite life time under a constant load in the spirit of the 
classical model of Coleman Q]. For natural fiber com- 
posites a so-called damage accumulation model has been 
developed, which simply assumes that the time deriva- 
tive of the accumulated damage depends exponentially 
on the external load history of the specimen S . 

In the present letter we study the creep rupture of fiber 
composites where the fibers have viscoelastic behaviour 



and the microscopic damage mechanism leading to creep 
rupture is the strain dependent breaking of fibers un- 
der the time evolution of the deformation of the system. 
Creep failure tests are usually performed under uniaxial 
tensile loading when the specimen is subjected either to a 
constant load a a or to an increasing load (ramp-loading) 
and the time evolution of the damage process is followed 
by recording the strain s of the specimen and the acous- 
tic signals emitted by microscopic failure events. In the 
present study we focus on the general aspects of creep 
rupture, i.e. the behaviour of the life time of the bun- 
dle as a function of the external load, its dependence 
on the range of load redistribution, furthermore, general 
aspects of the acoustic response of the bundle are consid- 
ered without fitting the theoretical results to any specific 
materials. 

In order to work out a theoretical description of creep 
failure of viscoelastic fiber composites we improve the 
classical fiber bundle model [Mpl which has proven very 
successful in the study of fracture of disordered mate- 
rials ]10|-pjj. Our model consists of N parallel fibers 
having viscoelastic constitutive behaviour. For simplic- 
ity, the pure viscoelastic behaviour of fibers is modeled 
by a Kelvin- Voigt element which consists of a spring and 
a dashpot in parallel and results in the constitutive equa- 
tion 



(T = j3e + Ee, 



(1) 



where (3 denotes the damping coefficient, and E the 
Young modulus of fibers, respectively. Eq. (|l]) provides 
the time dependent deformation e(t) of a fiber at a fixed 
external load a Q 
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where e D denotes the initial strain at t = 0. It can be 
seen that e{t) converges to <j q /E for t — > oo, which im- 
plies that the asymptotic strain fulfills Hook's law. 

If no fiber failure occurs Eq. (||) would fully describe 
the time evolution of the system. Motivated by the 
experimental observations of the acoustic response |2~l| l 
of fiber composites during creep, we introduce a strain 
controlled failure criterion to incorporate damage in the 
model: a fiber fails during the time evolution of the sys- 
tem if its strain exceeds a damage threshold Sd, which 
is an independent identically distributed random vari- 
able of fibers with probability density p(£d) and cumu- 



lative distribution P(£d) 



p(x)dx. Similar strain 



controlled breaking was recently used in Ref. j20[. Due 
to the validity of Hook's law for the asymptotic strain 
values, the formulation of the failure criterion in terms 
of strain instead of stress implies that under a certain 
steady load the same amount of damage occurs as in the 
case of stress controlled failure, however, the breaking 
of fibers is not instantaneous but distributed over time. 
When a fiber fails its load has to be redistributed to the 
intact fibers. As the simplest approach, we assume global 
load sharing jlO 15 -li]], i.e. after a failure event the ex- 
cess load is equally distributed among the intact fibers, 
and hence, at a certain strain s the load on the surviv- 
ing fibers of the number N s (e) can be cast into the form 
o-(e) = <j N/N s (e) = cr /(l - P(e)). The time evolution 
of the system under a steady external load <j is finally 
described by the equation 
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where the viscoelastic behaviour of fibers is coupled with 
the failure of fibers in a global load sharing framework. 

For the behaviour of the solutions of Eq. (||) two dis- 
tinct regimes can be distinguished depending on the value 
of the external load a : When a is below a critical value 
a c Eq. (^) has a stationary solution e s , which can be 
obtained by setting i = 



a = Ee,[l - P(e,)]. 



(4) 



It means that until this equation can be solved for e s 
at a given external load er , the solution s(t) of Eq. (||) 
converges to e s when t — > oo, and no macroscopic failure 
occurs. However, when a a exceeds the critical value a c 
no stationary solution exists, furthermore, i remains al- 
ways positive, which implies that for a > er c the strain of 
the system e(t) monotonically increases until the system 
fails globally at a time t 1 . 

In the regime a < a c Eq. (Q) also provides the asymp- 
totic constitutive behaviour of the fiber bundle which can 
be measured by controlling the external load a a and let- 
ting the system relax to e s . It follows from the above 
argument that the critical value of the load a c is the 



static fracture strength of the bundle which can be de- 
termined from Eq. (^) as a c = Ee c [l — P(e c )], where e c is 
the solution of the equation da /de s \ £ = 0, as shown by 
Sornette fl5|| . Since a (e s ) has a maximum of the value 
cr c at e c , in the vicinity of £ c it can be approximated as 



(5) 



where the multiplication factor A depends on the proba- 
bility distribution P. A complete description of the sys- 
tem can be obtained by solving the differential equation 
Eq. (||) . After separation of variables the integral arises 
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where the integration constant C is determined by the 
initial condition e(t = 0) = 0. 

The creep rupture of the viscoelastic bundle can be in- 
terpreted so that for a a < er c the life time (or the time to 
failure) of the bundle is infinite tf — oo, while above the 
critical load a > a c global failure occurs at a finite time 
tf, which can be determined by evaluating the integral 
Eq. (H) over the whole domain of definition of P(e). From 
the theoretical and experimental point of view it is very 
important how tf depends on the external load above 
<7 C . When a a is in the vicinity of cr c , i.e. a — er c + Acr , 
where Aa « <r c , it can be expected that the curve of 
e(t) falls very close to e c for a very long time and the 
breaking of the system occurs suddenly. Hence, the total 
time to failure, i.e. the integral in Eq. (^), is dominated 
by the region close to e c when Act is small. Making use 
of the power series expansion Eq. (||) the integral in Eq. 
(pf) can be rewritten as 
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which has to be evaluated over a small e interval in the 
vicinity of e c . After performing the integration it follows 

t f w (ao-^)- 1 / 2 , for 



cr > cr c . 



(8) 



Thus, tf has a power law divergence at a c with a univer- 
sal exponent — h independent of the specific form of the 
disorder distribution p(e). 

For the purpose of explicit calculations we consid- 
ered the case of a uniform distribution of the damage 
thresholds between and a maximum value e m , thus, 
p(ed) — l/e m and P{ed) — £d/s m - The stationary so- 
lution, the critical load and the corresponding critical 
strain can be obtained as <r = Ee[l— s/£ m ], a c = Ee m /4:, 
£c — s m /2, respectively. Finally, the solution of the inte- 
gral Eq. (g) taking the initial condition also into account 
can be cast into the implicit form 
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for cr < <j c (below the critical point), and 
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for cr > cr c (above the critical point). The behaviour of 
this analytic solution is illustrated in Fig. [j] for several 
different values of a . 
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FIG. 1. The analytic solution e(t) given by Eqs. (p|^o|) for 
several values of a below and above <r c . The critical strain 
e c and the time to failure tf for one example are indicated. 



The time to failure tf can be determined by setting 
e — e m in Eq. (JIG) , which results in the form 



-1/2 



(11) 



in accordance with the above general arguments. 

A further important general property of e(t) that can 
be deduced from Eqs. is that at the time to failure 
tf the deformation rate de/dt diverges. For disorder dis- 
tributions P(e) defined in a finite interval the exponent 
is universal de/dt « (tf — t)^ 1 ^ 2 . 

In order to obtain information about the gradual 
breaking of fibers during the creep process, in the exper- 
iments the acoustic signal emitted by breaking events in 
a short time interval is investigated. In our fiber bundle 
model the number of fibers Nb(t) which have been bro- 
ken up to time t can be determined as Nb(t) = NP(e(t)), 
and hence, its derivative provides the quantity 



1 dN b dPde _ p(e)Ee 
N~df ~ ~dedt ~ 
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which is a measure for the acoustic response. The be- 
haviour of Eq. jl2|) for the uniform distribution is illus- 
trated in Fig. where it can be observed that the acous- 
tic activity, i.e. fiber breaking, practically disappears in 
the plateau region of e(t) (compare to Fig. |l|), however, 
it diverges at tf due to the diverging deformation rate. 
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FIG. 2. The analytic solution for the breaking rate Eq. jl^ ) 
for several values of a Q below and above o c 



Since during a creep test e(t) is monitored from which 
de/dt can be calculated, furthermore, dN^/dt is mea- 
sured by means of acoustic emission techniques, Eq. ( p^ ) 
makes possible to determine experimentally the distribu- 
tion of the failure thresholds p(sd)- 

To complement the predictions of the analytic ap- 
proach Monte Carlo simulations of the failure process 
have been performed using global load sharing (GLS) 
and local load sharing (LLS) for the stress redistribution. 
The GLS simulation of the creep failure process of a bun- 
dle of N fibers proceeds as follows: (i) random breaking 
thresholds £j, i — 1 , . . . , N were chosen according to a 
probability distribution p, then the thresholds were put 
into increasing order, (ii) Since the fibers break one-by- 
one, the actual load on the fibers after the failure of i 
fibers is = a a N/ (N — i) where i — 0, . . . , N — 1, and 
the time between the breaking of the ith and i+lth fibers 



reads as U 
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nally, the time as a function of e can be obtained as 

t( £ i) — X)j=o^i( £ i) f rom which £j(t) can be determined. 
The time to failure tf of a finite bundle is defined as the 
time of the failure of the last fiber. To test the validity 
of the power law behaviour of tf given by Eq. (||) simu- 
lations were performed with various distributions in the 
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framework of GLS. The results are presented in Fig. || 
where an excellent agreement of the simulations and the 
analytic results can be observed. The macroscopic strain 
of the system e(t) and the acoustic response obtained by 
simulations was also found to be in agreement with the 
analytic results. 

To clarify how the damage process and the behaviour 
oft/ is affected by the range of interaction among fibers, 
i.e. by the range of load sharing we performed simu- 
lations with LLS on a square lattice of 200 x 200 sites, 
redistributing the load of the failed fiber on its nearest 
neighbors. The critical load a c was first determined as 
the static fracture strength of a dry fiber bundle with LLS 
assuming perfectly elastic behaviour for the fibers. Com- 
paring the results of the LLS simulations to the global 
load sharing results it was observed that above a c the fail- 
ure of the viscoelastic bundle occurs much more abruptly 
than in the case of GLS. 
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FIG. 3. The behaviour of the time to failure tf for uniform 
and Weibull distributions with two different Weibull moduli 
for the GLS case. All the three curves are parallel to each 
other on a double logarithmic plot with an exponent close to 
0.5 in agreement with the general result Eq. (ft). 



Summarizing, we studied the creep rupture of vis- 
coelastic fiber composites by enhancing the classical fiber 
bundle model. Varying the external load, two regimes of 
the creep process were revealed characterized by an infi- 
nite life time below, and by a finite one above the crit- 
ical load. The transition between the two regimes was 
found to be continuous for global load sharing, while it 
is abrupt for localized load redistribution, analogous to 
a second and first order phase transition, respectively. 
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FIG. 4. Comparison of 1/tf for the LLS and GLS cases as 
a function of the external load. 
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Varying a as a control parameter the two regimes of 
the creep rupture process are characterized by an infinite 
life time below a c and by a finite one above the critical 
point. The nature of the transition between the regimes 
in the global and local load sharing models can be char- 
acterized by studying 1/tf as a function of the control 
parameter a. In Fig. Hit can be observed that below the 
critical point, when no global failure occurs, 1/tf is zero, 
while above er c it takes a finite value for both LLS and 
GLS. However, the behaviour of 1/tf in the vicinity of a c 
is completely different in the two cases, for GLS the tran- 
sition is continuous, while for LLS 1/tf has a finite jump, 
analogously to a second and first order phase transition, 
respectively. 
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